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Abstract 

An overview of the computational results obtained from the NASA Langley developed 
unstructured grid, Reynolds-averaged Navier-Stokes flow solver USM3D, in support of the 
Ares-I project within the NASA’s Constellation program, are presented. The numerical data are 
obtained for representative flow conditions pertinent to the ascent phase of the trajectory at both 
wind tunnel and flight Reynolds number without including any propulsion effects. The USM3D 
flow solver has been designated to have the primary role within the Ares-I project in developing 
the computational aerodynamic data for the vehicle while other flow solvers, namely 
OVERFLOW and FUN3D, have supporting roles to provide complementary results for fewer 
cases as part of the verification process to ensure code-to-code solution consistency. Similarly, as 
part of the solution validation efforts, the predicted numerical results are correlated with the 
aerodynamic wind tunnel data that have been generated within the project in the past few years. 
Sample aerodynamic results and the processes established for the computational solution/data 
development for the evolving Ares-I design cycles are presented. 

Nomenclature 

AD AC Ares-I Design Analysis Cycle 

BDM Booster Deceleration Motors 

BTM Booster Tumble Motors 

Dia First stage diameter, ft 

Cl AD AC-2 A with full protuberances 

C2 ADAC-2A with simplified protuberances 

C4 clean ADAC-2A 

Caf axial force coefficient, Fx/q x S lc f 

Cn normal force coefficient, F z /q. x S lc f 

(Cn)c cumulative sectional normal force coefficient, F z /q. x S lc f 

C m pitching moment coefficient, pitching moment/q ro S re fX mrc 

Crm rolling moment coefficient, rolling moment/ q. x S lc fDia 

C p pressure coefficient, p-p.^/q^ 

CFD computational fluid dynamics 
CP center of pressure 
CPU Central Processing Unit 
DAC Design Analysis Cycle 
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FLT Flight 

F&M Force and Moment 

GN&C Guidance Navigation & Control 

LaRC Langley Research Center 

LAS Launch Abort System 

LH2 Liquid Hydrogen 

M Mach number 

NASA National Aeronautics and Space Administration 

NTF National Transonic Facility 

OML Outer Mold Line 

p local static pressure, psf 

p„ free-stream static pressure, psf 

PSWT PolySonic Wind Tunnel 

q x free-stream dynamic pressure, psf 

Re Reynolds number based on Dia 

RCS Reaction Control System 

RoCS Roll Control Systems 

S re f Reference area, ft 

SA Spalart-Allmaras 

SST Shear Stress Transport 

SRB Solid Rocket Booster 

UPWT Unitary Plan Wind Tunnel 

WT Wind Tunnel 

X/Dia Normalized axial distance measured from the main engine gimbal point 

XYZ Reference coordinate system 

X mrc Axial location of the moment reference center 

a angle of attack, degrees 

A incremental value 

O roll angle, degrees 


Introduction 

NASA’s Constellation program [1] was initiated to address the needs of U.S. access to space 
in the post space shuttle era with the mission to the moon, Mars and beyond. Unlike the space 
shuttle, where both crew and cargo are launched simultaneously, the architecture for the 
Constellation program proposed two separate launch vehicles: the Ares-I for the crew and the 
Ares V for the cargo. The Ares-I, the vehicle of interest in the present study, is a two-stage 
rocket that incorporates a modified five-segment solid rocket booster from the shuttle program 
while the upper stage is powered by the J-2X engine from the Apollo program. The major 
components of the vehicle and the proposed nominal flight trajectory are shown conceptually in 
Figure 1. The Ares-I major components consist of the Launch Abort System (LAS), Crew 
Exploration Vehicle (CEV) named Orion, the Upper Stage (US), and the First Stage (FS) Solid 
Rocket Booster (SRB) that includes the frustum, and the aft-skirt. 
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Launch Abort System 
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b)- Nominal flight trajectory [2], 

Figure 1. Ares-I major components and its nominal flight trajectory. 
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The Ares-I Design Analysis Cycle (AD AC) outer mold line (OML) definition has evolved in 
time, as shown in figure 2, with the progression of design refinements over the past few years. 
These design refinements are being incorporated to satisfy various discipline constraints such as 
the structural loads over the vehicle and on the individual protuberances, aerodynamic 
performance, flight guidance, navigation, and control (GN&C), etc. There are numerous 
protuberances on the vehicle, adding to the complexity of the geometry and the resulting flow 
physics of the numerical simulation, as well as the wind tunnel model design, fabrication and 
testing. For example, figure 3 shows a translucent front view of a typical Ares-I vehicle and a list 
of individual and the group of protuberances along with the abbreviations used in the present 
analysis. Such protuberances included the LAS nozzles, CEV umbilical cover, many others over 
the US and the FS such as the Roll Control System (RoCS), Reaction Control System (RCS), 
Liquid Hydrogen (LH2) feed-line, Systems Tunnel, Booster Deceleration Motors (BDM), 
Booster Tumbling Motors (BTM), hold-down wedges, and many other smaller ones such as 
cameras, antennas, and rings. 


ADAC1 


A 






i i 

Early 2006 


ADAC-2 (A101) ADAC-2B ( A 103 ) A105 



Early 2007 Early 2008 Late 2008 

Figure 2. Evolving Ares-I configurations. 
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2 Upper Stage Reaction Control Systems (2 US RCS) 

2 Upper Stage Roll Control Motors (2 RqQS) 

1 Umbilical Cover (1 UmkW 
4 Upper Stage UilMQM Motors (4 US Ullage) 

1 First Stage Tunnel (1 FSTUN) 

1 Upper Stage Tunnel (1 USTTUN) 

2 Booster Tumbling Motors (2 frustum BTM) 

2 Upper Stage Cameras (2 US CAM) 

1 Launch Abort System Tunnel (1 LAS tunnel) 

2 Upper Stage Ultra High Frequency Antennae (2 US UHF) 
8 Booster Deceleration Motors (8 BDMs) 

1 Upper Stage Liquid Hydrogen Fmdfm ((1 US LH2) 

4 Launch Abort System Nozzles (4 LAS Nozzles) 

1 Upper Stage Fill/Drain-Line (1 US F/D-line) 

4 Wedges 


1 US 

F/D-line v 


2 RoQS 


2 US RCS 



- lUmbil 


tunnel 



2 US 
UHF 


2 US 


1 LAS 


CAM 


4 Wedges 


Translucent Front View 

Figure 3. Typical Ares-I vehicle with protuberances. 


Complementary to ground based wind tunnel and flight tests, Computational Fluid Dynamic 
(CFD) methods are being extensively used to aid the design and the ascent aerodynamic database 
development of the proposed crew launch vehicle. One such CFD method is the NASA’s 
Langley Research Center (LaRC) developed Navier-Stokes flow solver USM3D [3] that has 
been designated to be the primary code to develop the necessary computational aerodynamic 
data within the Ares project. Two other widely used flow solvers, FUN3D [4] and OVERFLOW 
[5], have supporting roles to provide complementary results for fewer cases as part of a 
confidence building process to ensure code-to-code solution consistency. The initial data 
requirements from CFD was to provide the integrated longitudinal normal force distribution (i.e., 
sectional line-load) along the length of the vehicle for structural design and to provide an 
estimate of the Reynolds number (Re) effects increment from the sub-scale model testing in wind 
tunnels (WT) and full scale vehicle in flight (FLT). Such numerical data are very difficult to 
obtain experimentally, as it will be discussed later, at sub-scale model dimensions with all the 
protuberances attached. The nominal ascent flow conditions generally ranging over 0° < a < 8°, 
0.5 < M < 4.5, at various vehicle’s roll angles (O) over 360 degrees. Currently, there is no 
ground facility to support testing at flight Reynolds number for the Ares-I class of vehicles. 

At the inception of the Ares-I ascent aerodynamics project, an aerodynamic panel was 
formed [6] to provide the technical oversight over the data generation, results and analyses. The 
aerodynamic panel quickly became a venue for technical exchange of information between the 
experimental, CFD, and the management teams. This interaction fostered an expanded role for 
various technical groups and disciplines, and in particular the CFD team as shown in the diagram 
in figure 4. In addition to the initial two data requirements (i.e., sectional line-load and the WT- 
to-FLT Re incremental effects), the CFD data were used by various other elements of the project 
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such as the venting and debris analysis groups, and the load analysis teams for the Orion capsule 
and individual protuberances. Also, the availability of the CFD data prior to the start of wind 
tunnel testing was found useful in aiding the model design, protuberance fabrications, test matrix 
development, and instrumentation distribution in the regions of interest where major flow 
changes were occurring. Furthermore, the CFD data were utilized in a number of trade studies to 
address the first order aerodynamic effects from changes made to various protuberances such as 
geometrical modifications, repositioning, and/or the addition of a proposed new one to improve 
the vehicle’s performance. Examples of such trade studies included the height variations of the 
upper stage system tunnel and LH2 feed-line; and/or replacing of the upper stage pressure tunnel 
with the system tunnel, etc. Complementary to the experimental data, all the CFD results were 
also made available to the appropriate ascent aerodynamic database development team [7] as 
well as the uncertainty quantification group [8]. Although, the GN&C aerodynamic database was 
primarily developed from the wind tunnel data due to their consistency in acquiring large amount 
of data points in a timely manner, the available CFD results were used as a sanity check to 
anchor the measurements. In particular, the CFD results were found to be necessary for the 
extrapolation of the sub-scale wind tunnel data to flight Reynolds number, although such a 
practice has not been validated for the Ares-type vehicles. However, a limited set of high Re 
wind tunnel data, which reached to approximately one tenth of the appropriate flight level, were 
obtained from the Langley’s National Transonic Facility (NTF) at transonic speeds to aid 
anchoring the slope of the extrapolation curve at the lower end. Sample results that show a 
comparison between the computed and measured surface pressure coefficients for the ADAC-2B 
(A 103) at transonic speed for the high Re will be discussed later. 



★ Data difficult to obtain, if not practical, any other way than CFD. 

Figure 4. USM3D data development and applications in Ares-I project. 
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The objective of the present study is to provide an overview of the computational activities 
conducted at LaRC in support of the external flow aerodynamic analysis, without any propulsion 
effects, of the evolving the Ares-I vehicle at nominal ascent flow conditions up to the first stage 
separation (See figure lb). The emphasis will be placed on the correlations of CFD results with 
the available experimental data for selected configurations. This CFD paper complements two 
other aerodynamic overview reports with the focus on the Ares-I project from the aerodynamic 
panel perspective [6] and the experimental investigations [9]. Computational results are based on 
the three dimensional, Reynolds-averaged Navier-Stokes equations with the assumption that the 
flow is fully turbulent over the entire vehicle. In addition, efforts have been made to build upon 
these analyses by incorporating some power-on effects, in particular for the RoCS thrusters, 
using the same CFD methodology and code-to-code verification support from the OVERFLOW 
and Loci-Chem flow solvers [10 and 11]. Finally, the recent enhancements to the flow solver that 
facilitated such applications within Constellation program will be reported by Pandya et. al. [12], 

Computational Approach 

USM3D Flow Solver 


USM3D [3] is a tetrahedral cell-centered, finite volume, Euler and Navier-Stokes flow 
solver. Time integration follows the implicit point Gauss-Seidel algorithm, explicit Runge-Kutta, 
and local time stepping for convergence acceleration. The solver scheme allows various options 
for computing the inviscid flux quantities across each cell face. These include Roe’s flux- 
difference splitting (FDS), advection upstream splitting method (AUSM), flux vector splitting 
(FVS), and Harten, Lax, and van Leer with contact restoration (HLLC) scheme. MinMod and 
Superbee flux limiters have been incorporated in the flow solver to smooth out the flow 
discontinuities due to shock waves. Turbulence models include Spalart-Allmaras (SA) [13] with 
and without wall function, k-e [14], Shear Stress Transport (SST) [15], and algebraic stress 
models. In addition to standard slip/no-flow through (inviscid surface) and no-slip/no-flow 
through (viscous surface) boundary conditions, the method also includes a transpiration 
boundary condition that can be used to simulate propulsion effects. Although USM3D is used 
mainly for ideal gas simulation with constant y, options for variable y and limited real gas effects 
are also available in the code. 

Grid Generation 


The surface triangulations along with field tetrahedral volume grids were generated using the 
GRIDTOOL [16] and VGRID [17] software developed at LaRC. A rectangular box enclosing the 
vehicle is typically used to define the computational domain far-field boundaries. Each face of 
this rectangular box is located approximately several body lengths away from the configuration 
in the upstream, transverse, and downstream direction. The normal spacing of the first cell height 
above the vehicle surface is about 0.01 inches, in full-scale dimension, for the WT Re grid. Same 
surface grid definition is also used to generate the corresponding volume grid for the flight Re 
with a much smaller height for the first cell normal spacing to resolve the expected thinner 
boundary layer. As a result, the finer grid resolved boundary layer for the flight Re cases 
generally contained about 20% more volume cells than its counterpart grid designed for the 
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computations at wind tunnel Re. All computational grids used in the present analysis included a 
sting geometry representation of the WT model. 

Solution development, convergence, verification and validation 


Many efforts were made early on during the course of this activity to establish the 
appropriate parameter settings for the USM3D flow solver, such as the flux schemes, limiters 
and turbulence model to ensure the solution development remains stable over the nominal ascent 
range of the flow conditions [18, 19, 20 and 21]. This led to the development of an efficient 
script that would automatically generate appropriate input-stream parameters to the flow solver 
and submit series of cases such as a Mach or angle-of-attack sweep to a designated computer 
platform. Upon completion, a summary page pertinent to the solution convergence 
characteristics would be generated for quality assurance purposes. A typical summary page for 
the AD AC-3 (A 106) configuration is shown in figure 5. This summary page includes the 
designed statistical convergence data for the longitudinal and lateral aerodynamic force and 
moment characteristics as well as the flow solution residuals and the final surface flow property. 
It should be noted that all present computations have been performed at the NASA’s Advanced 
Supercomputer (NAS) complex located at the Ames Research Center. 
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Solution convergence is evaluated by monitoring the overall fluctuations in all six 
longitudinal and lateral aerodynamic coefficients, in particular the longitudinal aerodynamic 
coefficients (Cn, Caf, C m ) and the vehicle total rolling moment coefficient (Crm). A solution is 
considered converged when the fluctuations in the 4-component coefficients (Cn, Caf, C m , and 
Crm) become generally less than 0.5% of its average final value over the last 2000 iterations 
along with about 2-5 orders of magnitude drop in the solution residuals for density (i.e., 
Iog(r/r0)) and turbulent viscosity (i.e., Iog(tnu/tnu0)). In addition to this summary page, a series 
of post processing are performed on a given solution to generate data in various forms and 
formats useful for the project needs. 

A conventional approach, that has been refined over the course of this activity, is used for the 
verification and validation of the predicted numerical results at wind tunnel Re. This approach is 
shown in the diagram in figure 6. Within the Ares project, the verification of the computational 
results is viewed as the ability to demonstrate the correctness of the mathematical understanding 
of the underlying flow physics associated with such class of vehicles over the nominal ascent 
range of flow conditions. Whereas, the validation of the computational results is primarily based 
on how well the numerical predictions compared to the real-world experimental data. 



The verification process in the present activity is accomplished by comparing the numerical 
results that has been obtained from other Navier-Stokes based CFD flow solvers to ensure code- 
to-code solution consistency. Some representative comparison between the results from USM3D 
and OVERFLOW as well as FUN3D will be shown later. As part of this verification process, 
efforts were made to address the sensitivity of the computed results to various factors such as 
grid refmement/adaptation, turbulence models, flux schemes and limiters. Some representative 
results from such sensitivity studies will also be shown later. 

The validation process of the computational results is primarily made through comparisons 
with the available experimental data such as the surface Cp, the total force and moment 
coefficients, and the integrated longitudinal normal force coefficients. The validation effort was 
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conducted for two variations of the configuration; the clean vehicle and the vehicle configured 
with full protuberances. In addition to the total force and moment measurements, the clean 
configuration is primarily used to measure the surface Cp with sufficient density to facilitate 
accurate sectional line-load integration to validate the corresponding numerical predictions. A 
technique, developed by Samareh [22], was used to compute the sectional line-loads based on the 
circumferential integration of the surface pressure coefficients over small longitudinal slices of 
the vehicle. The computed sectional line-loads for the configuration with full protuberances can 
not be validated because it is difficult, or impractical, to instrument a sub-scale wind tunnel 
model of this size (i.e., typically a 1% to 2% of full scale vehicle) with a large number of 
pressure ports that are required for an accurate integration of the sectional line-loads. However, 
as it will be discussed later, limited surface Cp measurements have been obtained on the 
configuration with full protuberances at selected locations and are used towards a partial 
validation of the computational results but not for the line-loads. As a result, only the total force 
and moment coefficients measured on the model with full protuberances are used to validate the 
numerical predictions. Furthermore, a limited set of surface pressure data has also been obtained 
in the 0.5 < M < 1.1 range at high Re, for both the clean and the full protuberances 
configurations in the NTF for CFD solution validation purposes. Sample measured data at high 
Re and the correlation with the computational results will be discussed later. 

Results and Analyses 

ADAC1 


In the early 2006, the released OML definition of the ADAC1 configuration was used to 
extract the appropriate geometry description of the vehicle in IGES format. The IGES file was 
used to define the computational surface grid with different density, leading to the subsequent 
volume grid generation, and grid refinement study for the clean configuration. As part of the 
solution sensitivity to grid refinement study, four sets of grids with various densities were 
generated and the statistics, in terms of total number of cells in the computational domain and 
surface triangles, are shown in table 1. These grids were used to compute the flow at nominal 
ascent flow conditions using both the SA and the SST turbulence models. 


Table 1. Surface and volume grid statistics for the ADAC1 configuration. 


Coarse Grid (CG) 

Base Grid (BG) 

Fine Grid (FG) 

Super Fine Grid (SFG) 

Volume 

Cells 

Surface 

triangles 

Volume 

Cells 

Surface 

triangles 

Volume 

Cells 

Surface 

triangles 

Volume 

Cells 

Surface 

triangles 

~ 3xl0 6 

~ 28,000 

~ 5xl0 6 

~ 50,000 

~10xl0 6 

~ 90,000 

~19xl0 6 

~ 150,000 


The effects of turbulence model on computed Caf and Cn are shown in figure 7 for the base 
and super fine grids as a function of Mach number at a = 7°, WT Re. The figure also includes the 
corresponding measured data obtained from the Langley’s Unitary Plan Wind Tunnel (UPWT) 
[23] and the Boeing PolySonic Wind Tunnel (PSWT) [24], In general, the effect of grid 
refinement on the predicted Caf is shown to be small across the examined Mach range. Also, the 
results indicate that the effect of turbulence model on the computed Caf is relatively larger at 
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low supersonic and especially at low subsonic Mach number of 0.5. The effects of both 
turbulence model and the grid refinement appear to produce nearly a constant shift in the 
computed Cn for the entire Mach range, except at M = 0.5. With a given grid size, the SST 
model generally predict a higher Cn compared to the SA model. However, with a given 
turbulence model, the computed Cn decreases with grid refinement. In general, the measured 
data tend to fall within the scatter of the computational results associated with different 
turbulence models and various grid densities, and that the correlation between experimental and 
computational data is fairly good over the examined Mach number range. However, the solution 
development with the SA turbulence model was found to be generally more stable, robust and 
efficient relative to the SST model. Further discussion on the results and the effects from such 
studies of the grid refinement and turbulence models can be found in references 18 and 25. 




Mach Mach 

Figure 7. Computed effects from different turbulence model and grid refinement, and 
comparison with measured Caf and Cn for the clean ADAC1 configuration a = 7°, 

WT Re. 

Sample plots of the computed and measured sectional line-load [25] for two regions of the 
clean ADAC1 configuration are shown in figure 8, for M = 1.4, a = 7°, WT Re. The figure also 
shows the computed line-loads obtained from OVERFLOW and FUN3D that were generated as 
part of the code-to-code verification analysis. All computed results are based on the SA 
turbulence model and have been generated independently following a generally accepted 
approach for best practices associated with each method. In general, the computed line loads 
correlate reasonably well with one another as well as the measured WT data both in terms of 
overall trends and magnitudes. It should be noted that, with the proposed Ares-I trajectory, it is 
estimates that the maximum dynamic pressure would occur over the speed range of 1 .2 < M < 
1.8. As a result, special emphasis has always been placed in the aerodynamic performance 
analysis of various OML’s over this speed range and in particular at M = 1.6. 
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Figure 8.Computed and measured sectional line-load for the clean ADAC1 configuration, 

M = 1.4, a = 7°, WT Re. 


ADAC-2A (A101) 

In the early 2007, the OML definition for the AD AC -2 A (A101) configuration was released. 
A sample of the CFD results obtained at nominal ascent flow conditions for the ADAC-2A 
(A101) configuration [19] are presented in this section. Three geometric variations of the vehicle 
were modeled, as shown in figure 9, and they include the clean (ADAC-2A C, also referred to as 
C4), simplified protuberances (ADAC-2ASP, also referred to as C2) and the full protuberances 
(ADAC-2AFP, also referred to as Cl). Relative to the C4, the C2 configuration only 
incorporated the LAS nozzles and the umbilical cover geometries. Four sets of grids were 
generated for these geometries. The grid size statistics in terms of total number of cells in the 
computational domain and surface triangles are shown in table 2. The computational grids 
ranged from ~8xl0 6 cells for the C4 at WT Re to ~43xl0 6 cells for the Cl configuration at FLT 
Re. Computations were performed at both WT and FLT Re for Mach numbers ranging from 0.5 
to 4.0, at angles of attack of 0°, 4°, and 7°. In addition, the C2 configuration was used in a 
separate study to assess the viability of an existing grid adaptation technique for possible 
improvements in either the solution development or aerodynamic prediction accuracy. The 
results from this study led to no significant improvements, particularly for a milestone/time 
driven project environment, and thus were not pursued any further. 


Table 2. Surface and volume grid statistics for the ADAC-2A configuration. 


C4 (WT Re) 

C2 (WT Re) 

Cl (WT Re) 

Cl (FLT Re) 

Volume 

Cells 

Surface 

triangles 

Volume 

Cells 

Surface 

triangles 

Volume 

Cells 

Surface 

triangles 

Volume 

Cells 

Surface 

triangles 

~ 8xl0 6 

~ 220,000 

~ 21xl0 6 

~ 240,000 

~26xl0 6 

~ 350,000 

~44xl0 6 

~ 350,000 
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The computed and measured [23 and 24] Caf and Cn for the Cl configuration are shown in 
figure 10 as a function of Mach number at a = 7°, WT Re. The figure also presents the same 
computed coefficients obtained from OVERFLOW and FUN3D that were generated as part of 
the code-to-code verification process. The computed Caf values from all three flow-solvers agree 
reasonably well with one another, and slightly under-estimate the measured axial force 
coefficient particularly over the low to mid-supersonic speed range. The measured Cn values 
generally fall within the scatter of the computational results obtained from different flow solvers 
over the examined Mach range. The computed results for both coefficients clearly capture the 
overall characteristics of the measured data. 




j_i '‘I. iii .'.I ^ i_j 

12 3 4 5 

Mach 


Figure 10. Computed and measured Caf and Cn for the ADAC-2A (A101) with full 

protuberances (Cl), a = 7°, WT Re. 
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An effort was initiated during the course AD AC-2 A (A101) computational analysis in 
establishing the required procedure to isolate the contribution of each protuberance towards the 
vehicle’s total integrated aerodynamic forces and moments. A sample result from this activity 
that shows the Caf contribution, as a fraction of the vehicle’s total, from individual or a group of 
protuberances is presented in figure 1 1 for various Mach number at a =7°, and WT Re. For 
example, the computed results indicate that the majority of the Caf contributions from the 
protuberances are coming from the BDM’s, the stiffeners/kick rings, and the LAS nozzles. Also, 
the maximum Caf contribution from the BDM’s occur at M = 0.9 followed by a decreasing 
trends with increasing speed. Similar analysis has also been performed for all other computed 
coefficients. Such a data analysis was found useful in aiding the vehicle’s protuberances design 
from an aerodynamic performance perspective. In particular, the contribution of various 
protuberances to the vehicle’s rolling moment was found to be very important and a sample 
results will be shown later. 


Mach number variations, i.e., 
m050a7, i.,e., M=0.50, a=7° 



BDM 
I BTM 
RoCS 
USM 

I Nozzles 
Umbilical 
I RCS 
UHF 
USTTUN 
RCS 

F/D-line 

FSTUN 

rings 

Wedges 


Figure 11. Individual/group protuberance(s) contribution to the computed Caf for the 
AD AC -2 A (A101) Cl configuration, a = 7°, WT Re. 


The effects of protuberances and Reynolds number on the vehicle’s cumulative normal force 
coefficient (Cn)c, computed at M = 1.6, a = 7°, are shown in figure 12. The results indicate that 
the majority of the changes in the computed (Cn)c, due to the protuberances occur ahead of the 
frustum about mid upper stage at WT Re. In addition, the presence of the protuberances causes 
an increase in the vehicle’s overall normal force coefficients at these flow conditions. Similarly, 
the results on the right hand-side of the figure indicate a moderate drop in the computed (Cn)c 
with increasing Re from WT to FLT for the Cl configuration at M = 1.6, and a = 7°. Such a drop 
in the computed (Cn)c becomes noticeable at about mid upper stage and can be attributed to the 
expected thinner boundary layer at FLT Re. More detailed comparisons of USM3D with wind 
tunnel data and the results from other CFD codes for ADAC 2A configuration are discussed in 
Pao et. al. [25] with an emphasis on the importance of protuberances presence on the vehicle’s 
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rolling moment coefficient. It should be noted that this particular OML definition was frozen in 
preparation for the manufacturing time requirement for various components leading into the 
Ares-IX vehicle and the subsequent flight test. The computational activities were continued to 
support the pre- and post-flight test and the results were reported by Bauer et. al. [26]. 




a)- Effects of protuberances 


b)- Effects of Re 


Figure 12. Effects of protuberances and Re on Computed (Cn)c for the ADAC-2A 

(A101) configuration, M = 1.6, a = 7°. 


ADAC-2B (A1031 


This section presents a sample of the CFD results [20] obtained at the nominal ascent flow 
conditions for the ADAC-2B (A 103) configuration. This particular OML, released in early 2008, 
was subjected to an extensive wind tunnel testing [23, 24 and 27] that provided critical data to 
validate the computational results. The lessons learned, from the CFD analysis of the ADAC-2A 
configuration, on the importance of the protuberances and the roll angle on the vehicle’s overall 
rolling moment coefficient were also incorporated and expanded under this activity in support of 
the GN&C data requirement. As a result, the number of CFD solution development for this OML 
increased by nearly an order of magnitude to about 300 cases which resulted an increase in 
computer usage of about 1.5xl0 6 CPU hours. 

Both the clean and the full protuberances vehicle geometries were studied. The 
computational results and the correlations with the experimental data are presented here first for 
the clean configuration followed by the analysis for the vehicle with full protuberances. All 
computed results are based on the SA turbulence model unless noted otherwise. Grid statistics in 
terms of total number of cells in the computational domain and surface triangles are shown in 
table 3. 
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Table 3. Surface and volume grid statistics for the ADAC-2B configurations. 


Clean 

Full Protuberances 

Volume Cells 

Surface triangles 

Volume Cells 

Surface triangles 

- 20xl0 6 

- 270,000 

- 49xl0 6 

- 770,000 


At the outset, an extensive study was conducted on the clean configuration that addressed 
issues such as the solution sensitivity to grid refinement and the aerodynamic effects of modeling 
various height levels of small geometrical features such as the ramps, steps, and cavities. The 
results from the latter effort and the consultation with the experimental team identified the 
appropriate height levels for such small geometrical features from the 1%-scale model 
manufacturing standpoint as well as the resulting aerodynamic impacts. Also, the results from 
the grid refinement study revealed the adequacy of a base grid (-20x1 0 6 volume cells with 
~0.30xl 0 6 surface triangles) relative to a fine grid (-64x1 0 6 volume cells with -0.6x1 0 6 surface 
triangles) for resolving the major aerodynamic characteristics of the clean configuration. 

An efficient approach [28] was adopted for correlating the many measured surface pressure 
data points with the computational results and a sample comparison over the LAS/CEV and 
frustum regions are shown in figure 13 for M = 1.6, a= 7°, WT Re. In this approach the 
measured surface Cp data, shown in color-coded circles, are superimposed over the background 
computed-results that are plotted over the same contour level and range. The experimental data 
here are measured with 115 ports, that are distributed non-uniformly along the top model surface 
meridian, at a = 7°, M = 1.6, and for 15° roll angle increments. Thus any disagreements between 
the computed and measured data are highlighted by the changes in the color mapping. In general, 
the correlations between the two data are considered good. 



Figure 13. Computed and measured surface Cp on the clean ADAC-2B (A103) 

configuration, M = 1.6, a = 7°, WT Re. 


The effects of Mach number on the computed and measured sectional line-loads over the 
frustum region of the clean configuration is shown in figure 14 for a = 7°, O = 0°, WT Re. The 
results indicate significant changes occurring over the region due to compressibility effects. 
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Further analysis of the computed results revealed that such changes in the computed line-loads 
are primarily associated with a pocket of separated flow in the region and it’s characteristics 
changes with Mach number. The trends and the large variations in the measured line-loads due to 
compressibility effects have been predicted reasonably well by the method. 



Figure 14. Computed and measured line loads over the frustum for the clean ADAC-2B 

(A103), a = 7°, WT Re. 

An effort was also made to assess the effects of different turbulence models on the computed 
results for the clean configuration. In addition to the SA, the two-equation SST and the k-s 
turbulence models were selected for the study. The computations for this study were performed 
at a = 7°, and WT Re. The effects of different turbulence models on various flow and/or 
aerodynamic properties such as surface flow pattern and pressure distribution, longitudinal 
normal and axial force distributions, total skin friction and its components, etc. were analyzed 
[20]. Such analysis, although not discussed here, also incorporated the complementary results 
from other CFD codes as well as the available experimental data. As an example, the effects of 
turbulence models on the computed Caf and Cn for the clean configuration are shown in figure 
15 as a function of Mach number at a = 7°, WT Re. The variations of the computed Caf with 
different turbulence models were found to be about 10%, and that generally the SST model 
predictions were closer to those of the SA model. The results also indicated that the two-equation 
models predict higher Cn than the SA model at supersonic Mach numbers. In general, the results 
from this study, and those conducted earlier on the prior the OML configuration, revealed that 
the SA turbulence model was the most numerically robust, computationally efficient, and 
accurate in predicting the required aerodynamic data for this class of vehicles at the flow 
conditions of interest. 
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Figure 15. Computed Caf and Cn with different turbulence model for the clean AD AC- 

26 (A103), a = 7°, WT Re. 


The experience and the lessons learned from the grid generation on the prior configurations 
are used to develop a baseline grid for the ADAC-2B (A 103) with full protuberances. This 
baseline grid contained 49xl0 6 volume cells and 0.77xl0 6 surface triangles. Close-up views of 
the surface grid for three regions of the vehicle are shown in figure 16. The surface grids are 
generally clustered around the region where large changes occur in the local geometry such as 
the LAS/CEV juncture and/or any of the protuberance/body intersections. 



LH2 region 


CEV region 


Figure 16. ADAC-2B (A103) computational grid. 

The baseline grid was used for the computations at the nominal ascent flow conditions using 
the SA turbulence model. Unlike the axisymmetric surface geometry of the clean configuration, 
computations for the vehicle with full protuberances were performed at various roll angles, 
typically ranging from 0° to 360° with 30° increments. Sample results that show the effects of 
protuberances at d> = 0° and at other roll angles on the computed sectional line-loads at M = 1.6, 
a = 7°, WT Re are presented in figure 17 for the front portion of the vehicle. The corresponding 
geometry segment is included at the bottom to help with visual feature identification. The local 
effects on the computed line-loads due to protuberances and various roll angles are significant. 
However, such effects, as demonstrated by the umbilical cover (i.e., X/Dia ~ 21.6), appear to be 
relatively local and not longitudinally cumulative. As a result, the variations in computed line- 
loads diminish to an approximately the same level at about mid upper stage. 
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Figure 17. Computed effects of protuberances with roll angle on the ADAC-2B (A103), 

M = 1.6, a = 7°, WT Re. 

Limited surface Cp measurements have been obtained on the configuration with full 
protuberances that can be used for a partial validation the computational results. As discussed 
earlier, such limited sets of data are not dense enough over the vehicle to accommodate for an 
accurate integration of the sectional line-loads. These surface Cp data have been measured along 
the model top centerline, similar to the clean configuration, as well as some around a few major 
protuberances such as the LH2, the umbilical cover and the upper stage tunnel system. As a 
result, these surface Cp data along with the measured total forces and moments are the primary 
sources of the experimental data that have been used to validate the computational results for the 
configuration with full protuberances. 

Sample computed and measured surface Cp at a = 7°, WT Re, for the Mach numbers of 1.6 
and 0.9 for the configuration with full protuberances are shown in figures 18 and 19. The 
experimental data have been obtained from the PSWT [24], In addition to the top centerline, the 
correlations between the computed and measured data are also shown for two rows: one above 
and the other below the LH2. With exception of some differences over the frustum and the aft- 
skirt regions, the comparison of results indicates a reasonable agreement, both in terms of trends 
and magnitudes, between the measurements and the predictions at both Mach numbers. Although 
limited to a small region, such comparisons provide some degree of confidence in the numerical 
predictions of the surface flow characteristics. 
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Figure 18. Computed and measured surface Cp on the ADAC-2B (A103) configuration 

M = 1.6, a = 7°, WT Re. 



Figure 19. Computed and measured surface Cp on the ADAC-2B (A103) configuration 

M = 0.9, a = 7°, WT Re. 
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Prior to the recent Ares-IX flight test [26], there was no experimental data available for this 
class of vehicles at flight Re that could be used to validate the computational results. As a result, 
an effort was initiated to measure a limited set of high Re surface Cp, that reached to 
approximately one tenth of the appropriate flight level, from the Langley’s National Transonic 
Facility (NTF) at subsonic and transonic speeds (i.e., 0.5 < M < 1.1) to aid anchoring the slope of 
the extrapolation curve at the lower end. The NTF test was conducted using a 1%-scaled model 
for both the clean and with the full protuberances geometries at selected flow conditions of 
interest for various Reynolds number with the tunnel running on the standard air mode (i.e., low 
Re) as well as the cryogenic mode (i.e., High Re). The model was instrumented with pressure 
port distributions and flow conditions similar to the one tested in the PSWT that would provide 
an opportunity for a tunnel-to-tunnel data comparison at low Re. 


All CFD solutions were completed and the results were delivered to the project office prior to 
the start of the NTF test to aid the planning and test matrix development by the experimental 
team. The details of the test results from the NTF and the correlations with the delivered CFD 
data along with the measurements from the PSWT can be found in reference 27. In general, the 
comparisons of the computed and measured surface Cp, as well as the sectional line-loads were 
found to be in reasonable agreements for the clean ADAC-2B (A 103) configuration for both the 
low and high Re flow conditions. Similarly, the measured surface Cp distributions for the model 
with the full protuberances, although limited, were also found to correlate well with the 
numerical predictions. A sample comparison of the computed and measured surface Cp 
distributions for the top centerline along two rows adjacent to the LH2, for the high Re case at M 
= 0.9, a = 7°, is presented in figure 20. The results indicate good correlations between the 
computed and measured data in terms of trends and magnitudes for different regions with the 
exception of some disagreements over the centerline around the frustum and the aft-skirt region. 
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Figure 20. Computed and measured surface Cp on the ADAC-2B (A103) configuration 

at M = 0.9, a = 7°, high WT Re. 
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The computed and measured Caf and Cn for the ADAC-2B (A 103) configuration with the 
full protuberances are shown in figure 21 as a function of Mach number for various angles of 
attack. The wind tunnel measurements have been obtained from the PSWT for M<1.6 and the 
Langley’s UPWT for M>1.6. The increasing trends in the measured Caf, with increasing a, are 
predicted reasonably well across the examined Mach range, however, their magnitudes are 
slightly under-predicted. Similarly, the computed Cn over predicts the magnitudes of the 
measured data at all examined angles of attack, especially at low supersonic Mach numbers. 
However, the increasing trend in the measured Cn, with increasing angle of attack, is predicted 
reasonably well across the Mach range. 



Figure 21. Computed and measured Caf and Cn at various angles of attack for ADAC- 

2B (A103), WT Re. 


ADAC-3 (A 1061 

In early 2009, the latest A106 OML was released ahead of Critical Design Review (CDR). 
The most extensive computations so far in this project were performed on this particular OML to 
address an update for the aerodynamic effects due to various parameters such as angles-of-attack, 
Mach number, roll angle, at both the wind-tunnel and flight Reynolds number [21]. In response 
to the additional requirements and support of other elements within the project, there were about 
1,000 CFD solutions developed for this particular OML that consumed about 3x1 0 6 CPU hours. 
Sample results and trends in the computed aerodynamic data are presented here for the typical 
flow features, the sectional line-loads, the longitudinal aerodynamic coefficients, the rolling 
moment, and the center of pressure movement at representative ascent free-stream flow 
conditions and roll angle. Although no experimental data is used in this section, but based on the 
lessons learned from the prior verification/validation activities, it is believed that the present 
method’s viability for such application has already been demonstrated. 
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Table 4. Surface and volume grid statistics for the ADAC-3 (A106) configuration. 


ADAC-3 (WT Re) 

ADAC-3 (FLT Re) 

Volume Cells 

Surface triangles 

Volume Cells 

Surface triangles 

~ 70xl0 6 

~ l.OxlO 6 

~ 90xl0 6 

~ 1.0 xlO 6 


The early 2009 released OML for the A 106 configuration was used to generate the 
computational surface and the flow field volume grids. The complete configuration with full 
protuberances was modeled with approximately one million surface triangles and representative 
surface grid distributions are shown in figure 22 for three different regions. This surface grid was 
used to generate the corresponding volume grids of roughly 70x1 0 6 and 90x1 0 6 cells for the 
computations at the WT and FLT Re, respectively, as shown in table 4. The field grid used for 
the computations at flight Reynolds number incorporated a smaller height for the first cell 
normal spacing to resolve the expected thinner boundary layer and thus increased the total cell 
size by about 20x1 0 6 relative to the grid size used for the WT Re computations. 



LAS/CEV region LH2/Frustum region Aft-skirt region 

Figure 22. Computational surface grid for the ADAC-3 (A106). 


Sample low supersonic flow features computed at M = 1.6, O = 0°, WT Re, for a = 0°, and a 
= 8° are shown in figure 23. This Mach number is near the maximum dynamic pressure point 
along the ascent trajectory. The results show the computed surface pressure coefficient and the 
off surface Mach number contours. The off surface Mach contours, computed at a = 0°, show 
the expected initial bow shock forming ahead of the LAS nose followed by many other shock 
waves generated downstream over the body due to either the presence of a protuberance and/or 
surface geometry discontinuities. The low Mach number separated regions ahead and behind the 
ReCS and the RoCS seem to be small and do not appear to spill over onto the frustum. 

The computed flow features at a = 8° show a significant difference, as expected, between the 
top and the bottom half of the vehicle. The results show weakened shock waves on the windward 
side on the front LAS/CEV section whereas they become strengthened further down stream over 
the ReCS, RoCS and the frustum region. In an effort to highlight the separated regions, the 
portside view of the computed axial force component of the surface skin friction coefficients are 
shown in figure 24. In this figure, any reversed flow region is highlighted by the variations of the 

23 

American Institute of Aeronautics and Astronautics 


surface color contours. The results indicate a region of separated flow ahead of the LAS nozzles, 
over the CEY/ramp juncture and the first stage/frustum juncture, as well as around the BTM on 
the frustum. 
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Figure 23. Computed surface Cp and off-surface Mach contours for the AD AC -3 (A106) at 

M = 1.6, <E> = 0°, WT Re, portside view. 



Figure 24. Computed surface axial force component of the skin friction for the ADAC-3 
(A106) at M = 1.6, a = 8°, O = 0°, WT Re, portside view. 

The sectional line-load distributions computed at M = 1.6 and d> = 0° for various angles of 
attack are shown in figures 25. The results are presented with the focus on two sections of the 
vehicle: the LAS/CEV and the frustum section. Also, an approximate representation of the A106 
geometry profile is shown below the line-load plots to aid visual feature identification of axial 
locations. The results indicate, as expected, a nearly zero line-load distribution at a = 0°. At this 
angle-of-attack, any finite values of line-loads can be attributed to the geometrical asymmetry 
(i.e., protuberances) such as the LAS nozzles (X/Dia ~ 24). Note that the most significant load 
computed at a = 0° is occurring aft of the Interstage/frustum juncture where there are several 
asymmetrically positioned protuberances. As expected, the results generally show an increasing 

24 

American Institute of Aeronautics and Astronautics 





trend in the sectional line-load with increasing angle of attack with the exception of some 
variations over the frustum. 




a)- LAS/CEV region b)- Frustum region 

Figure 25. Effect of angle of attack on sectional line-loads for the AD AC -3 (A106) at M 

= 1.6, 0 = 0°, WT Re. 

The effects of roll angle on the sectional line-loads, computed at M = 1.6 and a = 8°, are 
presented in figure 26. The results at these selected roll angles indicate only a small change in the 
sectional line-loads over the LAS/CEV region, except around X/Dia ~ 19.8 where there are a few 
asymmetrically positioned protuberances. Also, the sharp spike in the sectional line-loads at 
X/Dia ~ 19.5 is due to the presence the first ring on the upper-stage. Unlike ADAC-2B (A 103), 
figure 17b, the smaller size umbilical cover for the present configuration eliminates nearly all the 
oscillations in the sectional line-loads (i.e., X/Dia ~ 21.6) with roll angle. The results also reveal 
a significant variation, due to roll angle, in the computed line-loads around the aft part of the 
upper stage and the frustum region. These large variations are caused by the flow interactions 
with many protuberances in the region such as the LH2, Upper Stage system and pressure 
tunnels, RoCS, and several other components. 
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a)- LAS/CEV region b)- Frustum region 

Figure 26. Effect of roll angle on sectional line-loads for the AD AC -3 (A106) at M = 1.6, 

a = 8°, WT Re. 


The effects of Mach number on the sectional line-loads, computed at a = 8°, O = 0°, are 
presented in Figure 27. The results show significant effects on the sectional line-loads due to 
Mach number. Such significant effects occur over the front part of the vehicle around the 
CEV/ramp juncture (X/Dia ~ 22.5), CEV shoulder (X/Dia ~ 21.7), and the first ring on the upper 
stage (X/Dia ~ 19.5). In general, the sectional line-loads show an increasing trend with Mach 
number over the CEV followed by a reversing trend after the shoulder (X/Dia ~ 21.8). Also, the 
maximum value of the sectional line-load occurs over the CEV (X/Dia ~ 22.4) at M = 4.5. The 
effects of the protuberances on the computed line-load over the frustum region appear to be most 
pronounced at M = 0.9 with a diminishing trend as the Mach number is increased. 


M=0.9, •1=0' 




a)- LAS/CEV region b)- Frustum region 

Figure 27. Effect of Mach number on sectional line-loads for the ADAC-3 (A106) at a = 

8°, WT Re. 
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The effects of vehicle’s roll angle on the computed Caf and Cn over the examined Mach 
range at a = 8°, WT Re, is shown in figure 28. The results indicate minimal roll angle effects on 
the computed Caf at all Mach numbers and the maximum level occurs at M ~ 1.2. However, the 
effects of roll angle on the computed Cn become more pronounced. The largest variations of the 
computed Cn with the roll angle occur at transonic and low supersonic speed range where it 
reaches to ~15%. Although not included here, further analysis of the computed flow features, 
sectional line loads as well as cumulative line loads indicated that the majority of the oscillations 
in Cn with roll angle at transonic and low supersonic speed range is due to the presence of many 
asymmetrically positioned protuberances over the aft-skirt. 
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Figure 28. Effect of roll angle on the computed Caf and Cn for the ADAC-3 (A106) at 

a = 8°, WT Re. 

A sample rolling moment coefficients computed at a = 8°, M = 1.6, and WT Re, for the 
ADAC-3 (A106) configuration is shown in figure 29. In addition to the total Crm for the vehicle 
(denoted as A 106), the figure also shows the corresponding contributions from the major 
protuberances such as al06_fp (i.e., A106 footprint: vehicle’s body with all protuberances 
removed), LH2, upper stage system tunnel (USTTUN), upper stage pressure tunnel (PressTun) 
Wedges, and etc. The results show that the maximum Crm for the vehicle occurs at d> = 0° and 
180°, and that the largest contributor to the total from the protuberances comes from the upper- 
stage LH2 tank fairing (its contribution is higher than the vehicle’s total). Other protuberances 
with a large contribution to the total Crm include the BDM’s, USTTUN, FSTTUN, RoCS, 
PressTun, ReCS, and al06_fp. Access to such data was found to be useful and was subsequently 
employed in several trade studies with the goal to reduce and thus improve the overall level of 
the vehicle’s rolling moment coefficients. Though not covered in this report, two approaches 
were investigated for such trade studies: one was through the circumferential redistribution of a 
selected protuberance (i.e., F/D-line, PressTun, BDM’s, etc.) by holding the axial location fixed 
and the other was to vary the protuberance’s height. Circumferential redistribution and the height 
variation of the selected protuberance were initially performed using the concept of superposition 
and the linear interpolation/extrapolation using the existing solutions. The accuracy of the results 
from such superposition and/or interpolation/extrapolation would later be verified for the most 
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promising options through a dedicated simulation study where the geometry changes were 
actually modeled to verify its viability. The results from such trade studies were distributed to the 
project team members for information towards a possible design changes or modifications. 
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Figure 29. Computed total Crm and the contributions from the major individual or 
group of protuberance(s) for the ADAC-3 (A106) at M = 1.6, a = 8°, WT Re. 

The ADAC-3 (A106) rolling moment coefficients and the center of pressure computed at a 
= 8° are shown in figures 30 as a function of the roll angle O, for various Mach number. The 
computed results indicate that the maximum Crm occurs at M = 1.6 for O = 0° and 180°. The 
results also show a vanishing Crm trend with increasing Mach number due to lower dynamic 
pressure at higher altitudes. The center of pressure (CP) is simply the ratio of the computed 
pitching moment to the normal force coefficients. It is a measure of the CP movement relative to 
the vehicle’s moment reference center (i.e., gimbal point of the main engine) in units of SRB 
diameter. The computed results show that there is a forward movement of the CP with increasing 
speed up to about Mach 3. The results also indicate that the variation of the CP movement in the 
axial direction, at transonic speeds (0.9 < M < 1.2), over the entire O-range is significant and that 
the level of CP-oscillations with roll angles is reduced with increasing speed to a negligible level 
at M = 4.5. It should be noted that sample computed force and moment coefficients for the 
ADAC-3 (A 106) configuration were subsequently correlated with the most recently obtained 
experimental data and the results are presented in reference 23. 
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Figure 30. Computed total Crm and CP at various Mach number for the AD AC -3 

(A106) at a = 8°, WT Re. 


The computed incremental effects due to Re on the sectional line-loads for various roll angles 
at M = 1.6, a = 8° are shown in figure 31. As expected, the results show that the most significant 
effects occur in the region where there is a large geometry change and thus the impact on local 
flow characteristics. Such regions, for the front portion of the vehicle, include the LAS nozzles 
(X/Dia —24.2 ), around the LAS ramp/CEV (22 < X/Dia < 23.2), around the beginning part of the 
upper stage system tunnel and some smaller protuberances such as the cameras, antennae, ring, 
and brackets (18.5 < X/Dia < 20). The incremental effects due to Re are more significant over 
the aft portion of the upper stage and frustum region where there are a number of geometrical 
changes. These geometrical changes include the frustum ramp, and many protuberances such as 
the LH2, F/D-line, four USM, two ReCS, two BTM, two rings, and the tail ends of the upper 
stage and pressure system tunnels. In general, the most significant incremental effects due to Re 
occur in the regions associated with the LH2 and F/D-line, the RoCS, and over the frustum ramp. 
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a)- LAS/CEV region b)- Frustum region 

Figure 31. Computed Re incremental effects on the line-loads at various roll angle for 

the ADAC-3 (A106) at M = 1.6, a = 8°. 
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The computed incremental effects due to Re on the sectional line-loads for various Mach 
numbers at a = 8°, d> = 0° are shown in Figure 32. Such incremental effects, though generally 
small, become noticeable over the LAS nozzles, around the LAS ramp/CEV juncture, aft portion 
of the upper stage, and the frustum region. Also note the presence of a significant effect around 
X/Dia ~21 at M = 0.9 that was found to be associated with a small pocket of separated flow 
caused by a shock/boundary-layer interaction. In general, the results indicate that these 
incremental effects due to Re diminish with increasing Mach number over the vehicle. 



a)- LAS/CEV region b)- Frustum region 

Figure 32. Computed Re incremental effects on line-loads at various Mach number for 

the AD AC-3 (A106) at a = 8°, O = 0°. 

The computed Caf and Cn for the AD AC-3 (A106) configuration at a = 8°, O = 0° for both 
the WT and FLT Re are shown in figure 33. It should be noted that these results are typical and 
that the general trends with respect to the Re effects remain the same for all examined roll 
angles. The results indicate a general reduction (~ 5% or less) of the computed Caf at flight Re 
across the examined Mach range. Similarly, the computed Cn at FLT Re also show a general 
drop in magnitudes, particularly over the low supersonic speed range (1.05 < M < 1.20) where 
they are reduced by as much as 7%. In an effort to find the cause of this reduction, the solution 
for M = 1.05 was selected for further analysis and will be discussed next. 



Figure 33. Computed Caf and Cn at WT and FLT Re for the ADAC-3 (A106) at a = 8°, 

0 = 0 °. 
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The computed flow features, consisting of surface Cp and off surface Mach contours, at a = 
8°, M = 1.05, O = 0°, WT Re and FLT Re, are shown in figure 34. Qualitatively, the effects of 
Re on the computed flow features appear to be very small over the front LAS/CEV region. 
However, such effects become significant over the aft upper stage and frustum region. For 
example, the computed surface Cp at FLT Re generally shows an expanded region on the surface 
that exhibit lower levels of Cp values (i.e., more negative) due to flow expansion in the aft region 
of the rings and the interstage/frustum juncture. In addition, the computed off surface Mach 
contours at FLT Re shows a trailing wake developing behind the camera on the leeward side of 
the frustum that doesn’t occur in the results obtained at WT Re. These effects can generally be 
attributed to the expected thinner boundary layer at FLT Re and the local flow interaction with 
any protuberances. 


Wind-Tunnel Re, M=1.05, a=8, <P=0° 



Figure 34. Computed surface Cp and off-surface Mach contours for the ADAC-3 (A106) 

at M = 1.05, a = 8°, <b = 0°, WT Re, portside view. 


Furthermore, the cumulative line-loads ((Cn) c) computed at a = 8°, M = 1.05, O = 0°, for 
both the WT Re and FLT Re, are shown in figure 35. The figure also shows the incremental 
effects on the computed (Cn) c (denoted by D(Cn) c) due to the Reynolds number. The results 
indicate a very small Re effects on computed (Cn) c up to about mid upper stage. In general, the 
computed D(Cn) c) results show that about one-half of the reduction in the computed Cn at FLT 
Re occurs over the aft portion of the upper stage and the other half takes place over the first 
stage. As a result, the changes in the computed line-loads over the frustum ramp, although 
locally significant, does not appear to contribute much to the overall reduction in the computed 
Cn at FLT Re. This overall reduction is found to be a cumulative effect due to the expected 
thinner boundary layer that results in a smaller viscous contribution to the normal force 
coefficient. 
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X/Dia 

Figure 35. Computed cumulative line-loads at WT and FLT Re and their increment for 

the AD AC-3 (A106), M = 1.05, a = 8°, O = 0°. 

A sample rolling moment coefficients computed at M = 1.6, a = 8° for both WT and FLT Re 
are shown in figure 36. With the exception at few roll angles, the computed Crm at FLT Re 
generally shows an increase relative to those computed at WT Re. This increase in the computed 
Crm at FLT Re becomes particularly important if it occurs at the maximum level and falls 
outside the roll control authority band. Such an increase in the computed Crm at FLT Re, though 
not validated, can have a significant effect on the design of an appropriate level of thrust 
requirements from the RoCS for the vehicle’s roll control during ascent. Such analyses have 
been undertaken within the project and the details are reported in references 1 1 and 12. 
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Figure 36. Effects of Reynolds number on computed Crm for the AD AC -3 (A106) at M 

= 1.6, a = 8°. 
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An effort was also made to evaluate the overall aerodynamic performance of the evolving 
Ares-I vehicle. Sample results computed at a = 7°, M = 1.6 and WT Re for Caf, and Cn are 
shown in figures 37, as a function of roll angle, for the selected configurations. With the 
exception of the computed Cn at O =180° and 210°, the results for the AD AC-3 (A 106) indicate 
a general increase in both Caf and Cn relative to the prior ADAC-2B (A 103) OML across the 
examined roll angles. Similarly the results computed at the same flow conditions for Crm, and 
CP are shown in figure 38. The results indicate that the AD AC-3 (A 106) configuration 
experiences an increased levels of maximum Crm at <3> = 0° and 180° relative to the prior design 
cycle OML. The variation of the computed CP with roll angle remains similar for both 
configurations, however, the overall magnitude has decreased for the AD AC-3 (A103) OML. 



Figure 37. Computed Caf and Cn for the ADAC-2B (A103) and ADAC-3 (A106) at a = 

7°, M = 1.6, WT Re. 



Figure 38. Computed Crm and CP for the ADAC-2B (A103) and ADAC-3 (A106) at a = 

7°, M = 1.6, WT Re. 

The notion of stake concept, as a passive flow control device, to reduce and therefore 
improve the rolling moment of the vehicle during the ascent phase of the flight trajectory was 
addressed. The effectiveness of the stake as a roll control device was initially established during 
the aerodynamic analysis of the ADAC-2B (A103) configuration and the computational results 
were subsequently validated against complementary data obtained from wind tunnel testing. In 
the initial study and the follow up validation effort, the proper axial and circumferential location, 
height, size, and even the dual use of the device were identified [21]. Although not discussed 
here, the lessons learned from this study were used to assess the effectiveness of the device for 
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improving the rolling moment of the AD AC-3 (A106) configuration. The configuration baseline 
grid was modified to incorporate the proposed strake geometry. This grid modification resulted 
in an increase of about 100,000 surface triangles and about 2xl0 6 volume cells over the baseline 
grid. The location of the proposed strake and a representative computational grid are shown in 
figure 39. The strake geometry can essentially be thought of as a cropped delta wing that its 
presence would induce a rolling moment effects which would counterbalance those of the LH2 
as a function of d> (see figure 29). The computed Crm, at a = 7° and 8°, M = 1.6 and WT Re, 
with and without the strake for the AD AC-3 (A106) configuration is shown in figure 40. The 
results with the strake indicate a significant reduction in the maximum Crm at both examined 
flow conditions. Representative comparisons of the computed results with the wind tunnel data 
for the configuration with and without the strake are reported in reference 23. Furthermore, the 
effectiveness of the strake as a roll control device for this configuration at flight Reynolds 
number was also verified. 



Figure 39. Shape, location and computational grid for the proposed strake on the 

ADAC-3 (A106). 


A106, no strake, a=7° 
A106, with strake, a=7° 
A106, no strake, a=8° 
A106, with strake, a=8° 


strake@270 



Figure 40. Computed Crm for the ADAC-3 (A106) with and without the strake at M : 

1.6, a =7° and 8°, WT Re. 
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Accomplishments 


An overview of the computational activities performed at the NASA Langley Research 
Center in support of the agency’s Constellation program and in particular the Ares-I project are 
presented. The numerical results are based on the turbulent Navier-Stokes formulation and have 
been obtained from the USM3D flow solver for the evolving design cycles of Ares-I 
configuration at representative flow conditions pertinent to the ascent phase of the trajectory. 
Computations have been conducted at both the wind tunnel and flight Reynolds numbers. A 
summary of the accomplishments based on the results and analyses are as follows: 

1- An efficient scripting process was developed for automating the solution generation and 
the associated convergence characteristic assessments to aid quality assurance, for a large 
number of configuration variations and flow conditions. 

2- Different forms of the geometry representations from a simple clean to the vehicle with 
the full protuberances have been modeled for various vehicles design cycles. 

3- Solution sensitivities to grid refinement have been addressed for a limited number of 
representative vehicle geometries. 

4- The effects of different turbulence models on the computed aerodynamic characteristics 
were investigated in detail for selected configurations and flow conditions. 

5- Computed flow solutions have been analyzed for aerodynamic features such as the 
surface Cp, sectional and cumulative line-loads, the overall longitudinal/lateral 
aerodynamic characteristic, center of pressure, and for their corresponding trends with 
respect to angle of attack, Mach number, Reynolds number and the vehicle’s roll angle. 

6- The aerodynamic performance increments from the sub-scale wind tunnel model to full- 
scale vehicle flight Reynolds number have been quantified for limited numbers of 
solutions obtained at wind tunnel and flight Reynolds number. 

7- An analysis process was established to isolate the aerodynamic contributions from 
individual or a group of protuberances towards the vehicle’s total. 

8- Demonstrated the feasibility of an approach used to estimate the aerodynamic impacts 
resulting from a possible repositioning/redistribution/addition/elimination of an 
individual or a group of protuberances on the vehicle’s overall aerodynamic performance. 

9- Established the viability of a new strake device to improve the vehicle’s rolling moment 
characteristics at both the wind tunnel and flight Reynolds number. 

10- As general practice within the project, all CFD solutions were completed and delivered in 
advance prior to the start of any experimental activities. Such practice was found to be 
useful for designing the wind-tunnel model hardware, its instrumentation sizing and 
strategic distribution, and test matrix development. 

11- Numerical predictions have been verified against the complementary results obtained 
independently from other flow solvers for solution consistency. 

12- Numerical predictions have been validated against the available experimental data that 
have been generated within the project over the past few years for various vehicles design 
cycles. 
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Concluding Remarks 


Feasibility of the present numerical method, as an engineering tool, has been 
demonstrated for the external flow aerodynamic analysis and design of such class of vehicles and 
flow conditions. The initial CFD data requirements, i.e. integrated line-loads and the 
aerodynamic effects resulting from sub-scale wind tunnel model to full-scale vehicle flight 
Reynolds number, were expanded during the course of this activity to include computational data 
analysis to support the need of other elements/groups within the project. These elements/groups 
include those supporting the debris and venting analysis, Orion loads, protuberance component 
loads, wind-tunnel model design/instrumentation, GN&C, and database development. In 
addition, USM3D data was used to study the impacts resulting from a possible repositioning, 
redistribution, resizing, addition, and/or elimination of an individual or a group of protuberances 
on the vehicle’s overall aerodynamic performance. For example, it was demonstrated that the 
addition of a new strake device would cause a significant reduction in the vehicle’s maximum 
rolling moment coefficient during ascent and its effectiveness was later validated by wind tunnel 
measurements. The USM3D results and analyses also indicated that the SA turbulence model 
was generally most accurate, robust, and numerically efficient as compared to other two equation 
models studied for the present application. 

Follow-on Efforts 

In spite of many efforts made during the course of these investigations to address the solution 
sensitivities to the grid refinement: a more systematic approach is presently underway to tackle 
the issue by providing a measure to assess the goodness (i.e., error/uncertainty band) of the A106 
baseline grid and the resulting predictions for the aerodynamic force/moment coefficients. This 
effort was made possible mainly due to the recent improvements and increase in the 
supercomputing power (Pleiades system) at the NAS facility located at the Ames Research 
Center. The grid densities used in this study ranged from 30x1 0 6 to 270x1 0 6 cells, and the 
computations are focused on conditions at flight Re with no propulsion effects. This study will 
explore the following issues: 

Establish a process to address the solution sensitivities to grid refinement for the Ares-I 
class of vehicles. 

Grid requirement for power off cases at flight Re. 

Asymptotic convergence behavior of the force/moment coefficients with grid refinement 
at various vehicle roll angles and Mach numbers. 
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